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^ . Abstract 



We consider four-dimensional de Sitter, flat and anti de Sitter 
branes embedded in a six-dimensional bulk spacetime whose dynam- 
ics is dictated by Lovelock theory. We find, applying a generalised 
version of Birkhoff 's theorem, that all possible maximally symmetric 
pq \ braneworld solutions are embedded in Wick-rotated black hole space- 

times of Lovelock theory. These are warped solitonic spaces, where the 
CN I horizons of the black hole geometries correspond to the possible posi- 

^+ ■ tions of codimension-2 branes. The horizon temperature is related via 

^^ I conical singularities to the tension or vacuum energy of the branes. We 

op ■ classify the braneworld solutions for certain combinations of bulk pa- 

rameters, according to their induced curvature, their vacuum energy 
and their effective compactness in the extra dimensions. The bulk 
Lovelock theory gives an induced gravity term on the brane, which, we 
argue, generates four-dimensional gravity up to some distance scale. 
d \ As a result, some simple solutions, such as the Lovelock corrected 

Schwarzschild black hole in six dimensions, are shown to give rise to 
self-accelerating braneworlds. We also find that several other solutions 
have self-tuning properties. Finally, we present regular gravitational 
instantons of Lovelock gravity and comment on their significance. 
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1 Introduction 

Combined cosniological and astrophysical data indicate that at least 96% of 
the actual matter content of the Universe has yet to be detected in particle 
accelerators. This is a correct statement provided we assume a homoge- 
neous universe described by Einstein's field equations. One fourth of this 
yet unseen component, dark matter, we can hope to discover in the Large 
Hardron Collider and is an expected, quite natural, cosmological signature 
of particle physics theories such as low energy supersymmetry (SUSY). For 
the remaining chunk, usually dubbed dark energy, that has cropped up more 
recently from supernovae data [1], the situation is a lot less clearer. Data 
is best "fitted" by assuming a tiny positive cosmological constant of magni- 
tude (10~^eV)'*. The mass scale associated to this energy is that of, roughly, 
what one expects for neutrino masses, a really tiny number in the cosmolog- 
ical arena. This cosmological constant is so small because the radius of the 
observed universe is by now huge, of size of the square inverse to the mass 
scale in question. It is dominant because dark matter is extremely diluted 
at cosmological scales of roughly lOGpc whereas super clusters of galaxies 
such as Virgo are at 20Mpc. 

Unlike dark matter, from the point of view of theoretical particle physics, 
the cosmological constant is rather problematic. For the "natural" value of 
the cosmological constant, from the point of view of Quantum Field Theory, 
is of the order of the ultraviolet cutoff we would impose for our quantum field 
theory, ranging from the Planck scale to the SUSY breaking scale, depend- 
ing on our theoretical prejudice. The fine tuning we impose to compensate 
this vacuum energy by some bare, gravitational cosmological constant is the 
biggest known discrepancy between theoretical expectations and our under- 
standing of experimental data. This is the cosmological constant problem [2] 
that was there already there [3] when the precision of existing data pointed 
towards a value, assumed for simplicity to be zero. 

By admitting the presence of a small, but not null, cosmological constant 
we add two additional problems to the "old" cosmological constant problem 
[3], namely, why is it so small and not zero, and why is its density now 
of the same order as dark matter density. Why therefore can we observe 
it now, i.e. the right time that it is possible for us to observe it in the 
first place? Dark energy models usually accompanied by scalar fields with 
suitable potentials do not solve this fine-tuning problem. They are also 
generically problematic due to enhanced radiative corrections associated to 
the scalar(s) [4] coupling to matter. Although general relativity is very well 
tested at solar system scales, both at weak and strong field regimes [5], it is 



not so at cosmological scales. In fact, to get an idea of the colossal difference 
in scales, take the quotient of the Hubble radius over the Astronomical unit 
(earth to sun distance) to get 10^^! Therefore, cosmological data, if we do 
not assume some extra source of unknown matter or cosmological constant, 
tell us precisely that gravity is modified at the infrared. 

Taking into account the above discussion, modifying gravity in the in- 
frared is a legitimate theoretical hypothesis that should be taken seriously. 
As yet, no convincing modification of general relativity has been found. 
Often the problem is experimental, constraints as for example for f{R) the- 
ories in the solar system and in galaxy clusters [6]. However, there are also 
theoretical difficulties since a large distance modification spoils our under- 
standing of asymptotic behaviour, which is very well-defined in General Rel- 
ativity (GR). Typically, such models [7,8] have ghost instabilities [9-11] and 
strong coupling problems [12] which are as yet not understood and therefore 
are no better than fitting data with a single, however small parameter, the 
cosmological constant. 

This of course does not mean there is no such consistent modifica- 
tion. Clearly one has to try more complex setups and some interesting 
toy models [13] have been proposed and are now being put to the test. By 
far, however, the most successful and popular modification is that of self- 
acceleration [14] present in the Dvali-Gabadagze-Porrati (DGP) model [8]. 
There, the small current acceleration of the Universe is understood as a ge- 
ometric effect originating from the (unfair) competition between the five- 
dimensional and four-dimensional curvature scales. Unfair, because the 
resulting length scale, the crossover scale, is of the order of the Hubble 
horizon today, Hq ~ 10~^^eU. At these length scales gravity becomes five- 
dimensional with a very low five-dimensional gravitational scale, roughly of 
lO^GeV and we enter a geometric acceleration phase. Unfortunately, this so- 
lution has a ghost in its linear perturbation spectrum [10,11]. Furthermore, 
the results are fuzzed by the presence of strong coupling problems [12] and 
either the theory loses linear predictability (which makes it useless cosmo- 
logically) or is completely unstable (see [15] for exact solutions and beyond 
linear order). To summarize, therefore, self-acceleration teaches us that 
we can in principle explain in a simple way the current acceleration of the 
Universe by a geometric effect of a codimension-1 braneworld, but the ge- 
ometry modification required to attain this result is not seemingly viable. 
Along these lines, recently [16] have given an idea on similar terms where 
the cosmological constant results from a higher order gravity correction in 
a codimension-1 braneworld. 

An even more ambitious proposal that emerged in braneworld models is 



that of the self-tuning and has to do with the "old" cosniological constant 
problem [3]. The idea was to find models where the tension of the brane, 
that is its vacuum energy, can be large without affecting the curvature of the 
brane and without fine-tuning of it with other brane or bulk parameters [17]. 
The best attempt to realize such a kind of model has been in six dimensions 
and in the case that the brane is of codimension-2 [18]. These conical branes 
have the special property that their tension merely induces a deficit angle in 
the surrounding bulk geometry and does not curve the brane world-volume. 
The brane curvature is induced from the bulk dynamics and is not directly 
related to the brane vacuum energy^. Several models with codimension-2 
branes were studied with various compactifications [19,20], however most of 
them had hidden fine-tunings or curvature singularities present. Although 
self-tuning vacua can be found, where indeed solutions of the same curvature 
correspond to a (continuous) range of brane vacuum energy, a cosmological 
constant problem resolution would require a dynamical selection mechanism 
of such vacua. This is because, as a rule, there exist nearby solutions of 
different curvature, which introduce the question whether the self-tuning 
solution is an attractor. 

In the present paper, we will examine a completely novel possibility of 
obtaining acceleration due to geometry as well as certain self-tuning prop- 
erties. The modified gravity theory that we will study is Lovelock the- 
ory [21] in six dimensions, which is the natural extension of GR in higher 
dimensions. The Lovelock theory in six dimensions has in addition to the 
Einstein-Hilbert term the Gauss-Bonnet combination. Although the lat- 
ter is a topological invariant in four dimensions, it becomes dynamical for 
higher dimensions and modifies the gravitational theory. The codimension- 
2 branes present in the generic vacua of this theory, can have interesting 
properties, similar to the ones stated in the previous paragraphs. We will 
find for example self-accelerating cases with non-compact internal spaces as 
well as self-tuning vacua for several compact and non-compact vacua. These 
examples open new possibilities for consistent self-acceleration and effective 
self-tuning which need to be considered in more detail in the future. 

In fact. Lovelock theory plays a crucial double role, not only in the bulk, 
but also for the brane junction conditions [22]. On the one hand, it provides 
geometric self-acceleration and novel solutions which are absent in Einstein 
theory, and on the other hand, it gives an induced gravity term on the 
brane [22,23], much like in DGP [8]. Here we must emphasise, that this 



^We emphasize here that we will not necessarily consider flat brane curvature, but also 
cases where a small brane curvature is allowed with a large vacuum energy on the brane. 



does not mean that we expect to have Einstein gravity on the codimension- 
2 braneworld, rather, it means that up to some scale, gravity "looks" four- 
dimensional as in the DGP model. Lovelock theory has the extraordinary 
geometric property to induce an Einstein-Hilbert term at the level of the 
junction conditions [23]. 

The structure of the paper is as followss: we will firstly present the 
general black hole solutions in six-dimensional Lovelock theory and then 
show how by double Wick-rotation we can obtain the most general axially 
symmetric braneworlds with maximally symmetric codimension-2 branes. 
Then we will scan through several classes of solutions and analyse the cases of 
interest, i.e. the self-accelerating and the self-tuning vacua introducing each 
time, as few "extra bulk parameters" as possible. We will briefly comment on 
the physical implications of some new regular instanton solutions of Lovelock 
gravity and we will finally conclude. 

2 Static black hole solutions 

Let us consider the six-dimensional dynamics of gravity with a bare cosmo- 
logical constant A and a Gauss-Bonnet (GB) term (see Appendix A for a 
more general system with in addition a gauge field coupled to gravity). The 
action of the system reads 



S= I <fx 
where 



{R + ccCgb) - 2A 



levrGe 



(1) 



-^GB = RmnkkR — '^RmnR + R , (2) 

is the Gauss-Bonnet Langrangian density, Gq the six-dimensional Newton's 
constant and a the Gauss-Bonnet coupling. 

It is straightforward to write down the Einstein equations of motion for 
the above action. They read 

Gmn - aHMN = — SvrGeAgfMAf , (3) 



2 

Bonnet contribution 



with Gmn = Rmn — 2^9 mn the Einstein tensor and the following Gauss- 

-^^GB Qmn — '^RRmn + ^RmkRn 
+4:RkmanR - "^RmkapRm ■ (4) 



Our ultimate goal is to find the maximally symmetric solutions of the 
above equations of motion. For this purpose, let us consider a four-dimensional 
space of maximal symmetry parametrised by k = 0, — 1, 1, with metric 

These spaces are four-dimensional flat space, hyperboloid or sphere respec- 
tively of curvature R[h] = 12k. A generalised version of Birkhoff's the- 
orem for (1) states, [24] (see [25] for the Lovelock version) that every six- 
dimensional spacetime solution of (3) having such four-dimensional maximal 
sub-spaces (5) is locally isometric to 



ds'^ = -V{r)dt^ + -— - + r^h^^dx'^dx" , (6) 



admitting therefore a locally timelike Killing vector (dt in the coordinates 
of (6)). The solution of the potential V for the equations of motion (3) is 
found to be [26, 27] 



l + eWl + 4a(a2-^ 



where the parameters appearing in the action (1) have been rescaled to. 



(7) 



a = 6a , levrGeA = 20a^ > (for dSa) , (8) 

a^ = -^2 < (for AdSa 



'6J 



The integration constant n is 

4^G6M 
^ = ^ 1 (9) 

where M is the AD or ADM mass of the solution and E^ is the area of 
the unit four-dimensional maximally symmetric subspace. Finally, e = ±1, 
giving rise to two distinct branches of solutions. The convention of the /x 
sign is chosen so that the gravitational mass is always /j, > 0. Indeed, as one 
can easily check by expanding the square root for large distances, the sign 
flip in front of ^ is necessary to match the Schwarzschild de Sitter solution 
behaviour for positive AD mass. The case where 1 + Aaa^ = is special, 
because the theory can be written in a Born-Infeld (BI) form [28]. 

Setting n = gives us the asymptotic vacua of the theory (1) which 
unlike Einstein theory are not unique (for given bare parameters in (1)). 



Then, we find that asymptotically we have an effective six-dimensional cos- 
mological constant 



levrGBA 
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•'■' a 



1 + e\/l + 4aa2 



(10) 



with a normalisation like in (1). We note in particular, that a gives an effec- 
tive cosmological constant in the bulk even without a cosmological constant 
term in the action, i.e. even with a^ = 0. The maximally symmetric space 
that the solutions asymptote to, depends on the sign of a and the branch 
of the solution, i.e. e. It is interesting to expand the above expression at 
a ^ in order to check if there is an Einstein theory limit for the above 
vacua or not 

levrGeAeff = (1 + e) - 20ea^ + 0(a) . (11) 

a 

Hence, we distinguish the following three cases 

• For e = +1, the solution, as can be seen by (11), does not have an 
Einstein theory limit as a ^ 0, although there may be a relevant dS^ 
or AdSQ Einstein solution mimicking (6). We name this branch the 
Gauss-Bonnet branch. The solution asymptotes AdS^ space for a > 
and to dS^ space for a < 0. Unlike what is argued in [26] this branch 
is not, at least, classically unstable for the reasons put forward by [29] 
(for a full study of this see [30] ) . 

• For e = — 1, the solution has an Einstein theory limit (as q ^ 0) as 
seen from (11) and therefore we call this branch the Einstein branch. 
The solution asymptotes AdSe space for a > 0, — l/(4a) < a^ < or 
a < 0, a^ < 0, to dSe space for a < 0, < a^ < l/(4|a|) or a > 0, 
a^ > and to Mq for a^ = 0. 

• For the special BI case l + 4aa^ = 0, as we can see from (10) we obtain 
an dS^ or AdS^ asymptotic solution which does not have an Einstein 
theory limit and no possible flat vacuum. However this is the only 
case that we have a unique vacuum. 

Let us now proceed into analysing the solution at hand (6) (we will 
mostly follow [31]). For the above metric solution (6), there are two possible 
singularities in the curvature tensor, the usual r = 0, and also a branch cut 
singularity at the (maximal) zero of the square root 



Whenever r^ is real and positive, this is the singular end of spacetinie (6). 
For the BI case 1 + 4aa^ = there is no such singularity. 

We will have a black hole solution if and only if there exists r = rh such 
that V{rfi) = and r^ > Tmax, where Vmax = max{0,rs}. Indeed the usual 
Kruskal extension 

{IT 

dv± = dt± -—— , (13) 

V [r) 

gives that (f+, r) and {v-,r) constitute a regular chart across the future and 
past horizons of (6). It is actually straightforward and rather useful to show 
the following: r = r^ is an horizon iff, 

e(2aK + r^)<0, (15) 

Pai^h) = with Pa{x) = —a X + Kx + aK X + eft . (16) 

In addition, one should make sure that the sign of V on each side of the 
roots of the above polynomial is such that (6) describes a black hole. The 
sign of Pa and V are not in general the same, depending on the signs of e 
and a. 

Therefore for e = +1 we have r^ < —2aK, i.e. , the event horizon is 
bounded from above. This means that for «; = there are no black hole 
solutions in this branch. Also note that Pa=oix) is just the usual polynomial 
for a Einstein black hole, and hence, a couples only to the horizon curvature 
K. Therefore, for e = — 1 and k = the horizon positions are the same as 
in the GR solutions. In fact, when e = — 1 we have similar structure and 
properties of the solutions (6) as their Einstein counterparts. 

3 Codiniension-2 braneworlds 

A standard procedure for Einstein theory to generate brane world solutions 
from black hole solutions [32, 33] is to perform a double Wick rotation for 
a black hole solution. The same will be valid for (1) and the black hole 
solution (6) at hand. We set t ^ iO and in addition make a further Wick 
rotation x^ — > it in the metric h^i, (5) so that it becomes of Lorentzian 
signature 

Then, the maximally symmetric spacetime sections correspond to four-dimensional 
Minkowski, AdS^ and d54 for k = 0,-1,1 respectively with curvature 



R[h] = 12k. The solutions are now of manifest axial symmetry with de 
as the angular Killing vector, 

ds^ = V(r)de^ + —— + r^K^dx'^dx'' . (18) 

V[r) 

The six-dimensional spacetime has the correct symmetries to describe a 
maximally symmetric four-dimensional brane world. The staticity theorem 
invoked for black hole spacetimes [24] tells us now that axial symmetry comes 
for free and need not be imposed for resolving the system of equations. The 
solutions (18) are therefore the general solutions describing maximally sym- 
metric branes, something that has been overlooked up to now even in the 
case of Einstein theory. It is interesting to note that in the case of Einstein 
theory in four dimensions, we obtain in this way the general maximally sym- 
metric gravitating cosmic strings, which can differ drastically from their flat 
counterparts (see in particular, the cosmic string solutions in AdS presented 
in [34]). 

The positions of the horizons r/j will be the endpoints of the internal 
space and the solutions will have meaning if we keep the spacelike regions 
of the previous black hole solutions, i.e. the ones with V{r) > 0. Let us 
suppose that we keep the space-like region between two horizons r_ and 
r+, with r_ < r+ [the subsequent discussion applies of course also in the 
case that we have only one horizon and we keep the side which is spacelike]. 
At these endpoints of spacetime, which are also the fixed points of the axial 
symmetry, one can in general put branes of dimension four, in other words 3- 
branes. Since the brane solutions we have found are maximally symmetric, 
the brane can only carry some two-dimensional Dirac charge, i.e. pure 
tension. Taking x^ = const, and expanding around the zeros of V we get. 



ds'^[iv;i)pide' + dpi, (19) 

with the Gaussian Normal radial coordinate 



P± = J^y, , (20) 



'r± 



which is well defined in all cases with V' ^ 0. The case of V' = 
needs special attention as we will see later (Sec. 4. 1.3 and Sec. 4. 2. 3). In 
this coordinate system, the brane energy momentum tensor is T^™"-^ = 

Sfiu^^^' = Sf^i, 2^"^ ; with S'^ = —T±6'i^, where T± are the brane tensions. If 



the angular coordinate has periodicity 9 G [0, 27rc), then the deficit angles 
which are induced at the two brane positions are 5± = 27r(l — f3±) with 

/?± = ^l^±|c . (21) 

The angular periodicity c is an arbitrary topological integration constant of 
the solution and can be varied to generate physically different brane world 
solutions. Let us note here, that in a compact model, we can always use 
the freedom of choosing c to set the deficit angles of one of the branes to 
zero, e.g. (3- = 1, thus, obtaining a "teardrop" compact space with only 
one brane of /3+ / 1. 

From the Einstein equations (3) supplemented by the brane tension 
terms, one can separate the distributional Dirac parts and write down in- 
duced Einstein equations for the branes. These brane junction conditions 
are [22,23] 

2vr(l - /3±) (-7^, + AaG%'^) = 87rGeS^, , (22) 

where 'yf^,^ = r'^h^^, is the induced metric on the branes with curvature 
R[j±] = l2K/rl = 12kHI, and G'^^^ = -SkHIj,,^ is the induced Einstein 
tensor. It is important to emphasize here that r^. depends on the geometric 
parameters of the bulk solution (18) namely the mass n, and the bare pa- 
rameters a and a^. It is also effectively the warp factor of the brane. The 
induced Newton's constant on the two branes can be determined from (22) 
to be 

^' = 4vra(l-/3±) ' ^^^^ 

Note that in order to have positive induced Newton's constant, we should 
have have angle deficit {P± < 1) for a > and angle excess (/3± > 1) for 
a < 0. The two dimensional warped space in the coordinates 9 and r has in 
the latter case the shape of a pumpkin whereas in the former case the shape 
of a football (rugby ball to be more precise). The important information 
coming from (22) is that the parameters a and (3 invoke the relation and 
the possible hierarchy between the bulk and four-dimensional scales. 

Substituting the G!^*^ back in (22) we find a relation between the Hubble 
parameters on the branes and the action parameters 

2.a_,„(-L,««|).l^r^ ,24) 

which can be further simplified, if we substitute the deficit angle from (23) 



and solve for Ha. 



,2 _ 1 87rG± ± 



KHi = + ^r± . (25) 

^ 2a 3 ^ ' 

The above equation is very important and relates the curvature on the brane 
H^ to its sources, namely the brane tension and the Gauss-Bonnet coupling. 
We see that the junction conditions tell us that the effective expansion H± 
is in one part due to the Gauss-Bonnet induced cosmological term and in 
another part due to the vacuum energy of the brane. Since H± = l/r±, we 
remind the reader that the Hubble parameter can be expressed as a function 
of parameters appearing in the soliton potential (7). Then, according to 
the specific explicit solution, the Hubble parameter is constrained and the 
above relation can give bounds on the two previously mentioned sources of 
curvature. 

A final comment before discussing the various brane world solutions has 
to do with the induced effective Newton's constant G^ appearing in (23). In 
a usual Kaluza-Klein dimensional reduction, the effective Newton's constant 
is obtained after substitution of the graviton's zero mode wavefunction in 
the action and the integration of the extra dimensions. This integration, 
actually, defines a relation between the effective Newton's constant on the 
brane, the higher dimensional Newton's constant and the volume of the 
internal space (see Appendix C for the volume calculation for the general 
brane world models that we will discuss). However, because of the presence 
of the branes and the Gauss-Bonnet bulk dynamics, the graviton dynamics, 
as perceived on the brane, have peculiarities and four dimensional dynamics 
with an induced Newton's constant G^ can be operative, even for internal 
spaces of infinite volume [8]. This is the interpretation one should give to 
the effective G^ appearing in (23) which can be different from the Newton's 
constant obtained via the volume calculation in Appendix C. For a finite 
volume element the actual Newton's constant perceived by the brane ob- 
server depends on the scale on which we probe gravity on the brane. This 
is analogous to what happens for the case of the DGP model [8] when we 
embed it in a Randall-Sundrum setup. 

4 Braneworld solutions 

In this section, we will analyze some particular cases of brane-world solu- 
tions, keeping only some of the parameters in the potential (7) each time, 
which give some interesting brane world examples. In each subcase, we will 
present the important features introducing the fewest parameters possible. 
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Figure 1: The horizon positions r_ and r+ as a function of the black hole 
mass /i for the Gauss-Bonnet branch e = +1 and for a < 0, where the vacua 
are ^5*4. For /j, = /j^c we obtain the degenerate horizon case. 

What will interest us in particular are the zeros of V which will correspond 
to possible brane locations. We will continue to refer to roots of V as hori- 
zons although this term is strictly speaking correct only for the black hole 
solutions (6). 

4.1 Zero bare cosmological constant 

Let us first choose zero bare cosmological constant a^ = 0. Then in all cases 
of this class, the horizon position is given by the solutions of the algebraic 
equation 



x-^ + Ax + M = . 



(26) 



The roots of the above equation are analyzed in Appendix B. In the present 
section, we will in addition apply the constraints (14)- (16) and the require- 
ment to be in a spacelike region with F > 0, in order to find all the possible 
vacua. 

Firstly, it is evident that for a^ = there are no k = flat vacua since 
the polynomial (16) has no roots^. The other cases depend on the branch 
and the four-dimensional curvature k. 

4.1.1 The dSi vacua 

Vacua with k = +1, which are dS^, exist both in the Gauss-Bonnet and 
the Einstein branches. As we will see, the Gauss-Bonnet branch vacua are 



■^This is true for solutions of the specific ansatz (18), which does not include the un- 
warped flat case. As we will see in Sec. 4. 1.5, one unwarped flat vacuum for a^ = exists. 
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Figure 2: The horizon position r^ as a function of the black hole mass ^ for 
the Einstein branch e = — 1, where the vacua are dS^. The point fi = fig 
corresponds to a singularity r^ and is excluded. 

compact with respect to the (r, 6) sections, while the Einstein branch vacua 
are non-compact. 

For e = +1, the only solution satisfying all the above requirements is 
for a < and < fi < fic = 2|ap'^/(3\/3). Then we have a double root 
structure, which means that the corresponding brane world solution is of 
compact (r, 6) sections. It is important to note that in this case the horizon 

< v^|a|/3 and i/|a|/3 < r+ < '~~ 



\a\ 



positions are bounded from above r_ 
so the corresponding Hubble parameters H± = l/r± are going to be bounded 
from below. When in particular n = Hc, the two roots become degenerate 
and y/^ = 0. The solution deserves special attention and will be discussed 
later on. The plot of the horizon (or brane) positions as a function of // 
is given in Fig.l. It is important to stress that the relevant black hole 
solution corresponding to this soliton, has the characteristics, of a de-Sitter- 
Schwarzschild like black hole studied in [31]. Note that the bigger the a 
the smaller the effective cosmological constant. Here we have a pure Gauss- 
Bonnet soliton (or black hole) in the sense that the a — > limit is singular. 
For e = —1, i.e. in the Einstein branch, we have two classes of solutions. 
The first class is for a > and /i > 0. The single horizon of that case 
can take all positive values. The second class is for a < and ji > ^^ = 
\/2|an^- In this case the horizon position is bounded from below as r^ > 
^j2\a\. Furthermore, in the latter case the singularity (12) exists and is 
hidden behind the horizon. In the limit /i — > /x^, we have r^ -^ Tg and the 
brane position becomes singular. We should therefore exclude the fi = fj,s 



point from the physical parameter space. In both cases, the brane Hubble 
parameter can be as small as desired. The plots of the horizon position 
as a function of fi in the above cases is given in Fig. 2. It is interesting 
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Figure 3: The horizon position rh as a function of the black hole mass fi for 
the Gauss-Bonnet branch e = +1 and for a > 0, where the vacua are AdS^. 

to note that both cases here correspond to a Gauss-Bonnet corrected six 
dimensional Schwarzschild black hole solution with a > or a < 0. 

4.1.2 The AdSi vacua 

Vacua with «; = —!, which are ^^5*4, exist only in the Gauss-Bonnet branch. 
In the Einstein branch, for all the cases satisfying the constraint (15), the 
region shielded by the horizon is timelike and therefore not suitable for our 
compactification. 

In the Gauss-Bonnet branch e = +1, we have only non-compact brane 
world models, all of them for a > 0. The solutions satisfy the criteria 
(14)-(16) for -/Xc < /^ < Aib = 5\/6a3/V9- The singularity (12) exists for 
< IJ, < fj,b but is always hidden behind the horizon. The plots of the 
horizon position as a function of /n in the above cases is given in Fig. 3. 

4.1.3 The Nariai dS^ vacuum 

The degenerate case for k = +1 in the Gauss-Bonnet branch e = +1 deserves 
special consideration. This happens when the mass parameter reaches its 
maximal value // — > /Xc and so the two horizons become one, r_ — > r+ — > 
y|a|/3. The internal space then appears to collapse, however one should 
look as for the Nariai metric to a different coordinate system [39]. So, setting 
^ = r_ /r+ , we make the coordinate transformation that consists of an affine 
transformation 



r = -f[{l-OR + l+C] 



2r+ 
(1-0 



e 



(27) 
(28) 
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The degenerate case corresponds to the hmit ^ ^ 1 and as can be seen from 
(18), it is the limit that the four-dimensional part of the metric is rendered 
unwarped. In these coordinates R £ [—1, 1] and Q G [0, 27rC) with 

+ -C . (29) 



(1-0 

The affine transformation blows up a point of the coordinate system (r, 6) to 

a well behaved internal space in the {R, 0) coordinates. Let us now define 

the modified potential 

AV 

with the parameters fi and |a| as functions of r-|_ and ^ given by 



/ = T^^^2 ' (30) 



Then, substituting everything in the metric (18) we obtain the blown up 
metric 

ds^=rl(fde'^ + '^ + ^h^^dx^'dx''^ . (32) 



The Nariai limit ^ ^ 1 gives a non-singular limit for the modified potential 
5^ 



/ -^ 1(1 — R?). Therefore, we obtain a non- vanishing internal space 



d.2 = ( M 



2 



^(1 - R^)d<d^ + 3-^^ + h^udx^dx' 
5 1(1-^2) 



(33) 



Note that, after this coordinate transformation, the deficit angle from 
(21) can be expressed as a function of the modified potential as 

l5+=P^ = \\dRf\C = \c . (34) 

2 5 

Thus, in this limit the deficit angle is independent of a and depends only 
on the periodicity C of the angular coordinate. 

4.1.4 Zero black hole mass fi = 

Let us now suppose that also the black hole mass vanishes /j, = 0. Then, the 
potential has the simple cosmological form 

2 

V{r) = K+^{l + e), (35) 
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from where we can see that the Einstein branch e = — 1 can only give the 
trivial flat six-dimensional solution, for k = +1. The Gauss-Bonnet branch 
e = +1 on the other hand has non-trivial solutions. The flat k = case 
is excluded since it has no horizon, but there are AdS4, and dS/^ solutions. 
Firstly, the k = — 1 case can have brane world solutions as long as a > 0. 
Then the spacelike region that we use for the internal space is non-compact. 
This case corresponds to the point of Fig. 3 where the curve intersects the 
/i = axis. 

On the other hand, for the k = +1 vacua we can have an horizon for 
a < and obtain a compact brane world model which is part of a six 
dimensional de-Sitter space. This case corresponds to the point of Fig.l 
where the upper curve of r+ intersects the ^u = axis. Furthermore, since 
the space has no singularity at r = 0, we can extend the radial coordinate 
to r < and consider the region of — ylaj < r < yr^- ^^^^ internal space 
is symmetric around r = 0, thus we have Z2 symmetry around the equator 
of the internal space. The reintroduction of a black hole mass /i 7^ breaks 
this symmetry since it introduces an r = singularity. 

4.1.5 Flat vacuum 

As noted before, the black hole ansatz (18) does not admit flat vacua for 
a? = 0. However, this ansatz, as we noticed also in the Nariai vacuum, 
has limitations, since it describes only warped four-dimensional metrics. 
Looking for unwarped solutions, it is easy to see that the following flat 
vacuum is a trivial solution of the equations of motion (3) 

ds^ = ri^^.dx^dx'' + dp^ + p^dO"^ , (36) 

where 6 G [0, 27rc). It is non-compact and exists for any sign of a. The 
equations (23), (24), (25) continue to hold for k = 0. To compactify such 
flat vacua, a gauge field flux has to be added along the lines of Appendix A. 

4.2 Non-zero bulk cosmological constant 

Let us now switch on the bulk cosmological constant, i.e. a? 7^ 0. The 
polynomial (16) is more complicated to solve, therefore we will study some 
special cases and focus on as few parameters as possible. Namely, the case 
with zero black hole mass ^ = and also the point l+4aa^ = in parameter 
space where the theory can be written in a BI form. Finally, we will provide 
all the flat (k = 0) vacua. 
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-1 





l/(4|a|) 







Figure 4: The horizon position Vh as a function of the bulk cosmological 
constant a^ for the dS^ vacua with // = 0. The solution on the left is for 
a < 0, where the Einstein and Gauss Bonnet branches meet at the BI point, 
represented by dot. The solution on the right is for a > 0. 



4.2.1 Zero black hole mass /x = 



If the black hole mass vanishes /i 
logical form 



0, the potential has the simple cosmo- 



V{r) = K + — [1 + e\/l + 4aa21 
2a 



(37) 



It is obvious that the only brane world solutions that we can obtain in the 
present case are the curved ones k = ±1. 

Let us first discuss the dS^ {k = +1) solutions. For a > 0, we see from 
(15) that only the Einstein branch e = — 1 has a solution with a horizon. 
Furthermore, in this case, in order to have r^ > 0, the bulk cosmological 
constant should be positive, i.e. a^ > 0. For a < we can have solutions 
both in the Einstein and the Gauss-Bonnet branch. In the Einstein branch 
e = — 1 the horizon is r^ > for positive bulk cosmological constant with 
< o^ < l/(4|a|). On the other hand, in the Gauss-Bonnet branch e = +1 
there are solutions for any a^ < l/(4|a|). At the BI limit a^ = l/(4|a|) 
the two branches merge to a unique dS^ solution. In all the above cases 
the spacelike region is for r < rh and the internal spaces are compact. As 
in Sec. 4. 1.4, we can extend the radial coordinate to r < and consider the 
region of —rh < r < r^, since the space has no singularity at r = 0. The 
plots of the horizon position as a function of /x in the above cases are given 
in Fig.4. 

The AdS^ (k = —1) solutions are the same as above if we substitute 
a^ -^ —c?. Thus, the curves in Fig.4 should be the mirrors of the k = \ case 
with respect to the a^ = axis. In addition, the spacelike region for these 
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Figure 5: The horizon position r^ as a function of the black hole mass ^ for 
the Born-Infeld case E = —1 and for a > 0. The left graph is for dS^ vacua 
and the right for AdS4 vacua. 

solutions is for r > rh, thus all the AdS^ vacua are for non-compact internal 
spaces. 

4.2.2 The Born-Infeld vacua 



For the dS4 and AdS4 vacua, it is easy to see that V is monotonic for E = —1 
and has an extremum at Te^tr = ilctlfJ-/^)^ if -E' = +1- 

The vacua obtained for 1 + Aaa^ = deserve special attention. As 
discussed in the Sec. 2, these vacua do not have an Einstein limit. In a 
certain way they correspond to the strongly coupled limit of the Gauss- 
Bonnet term in the action (1), since at the linearized level the combination 
(1 -|-4aa^) multiplies the perturbation operator. Furthermore, as we saw in 
the previous section, for these vacua the branches for e = +l,e = — 1 merge. 
Let us now switch on a mass parameter ^u > in the potential and write it 
as 



Vir) 



r^ M 

K+ — + ^ 

2a Jr 



(38) 



where the integration constant M = ^q- — replaces ^. The potential 
(38) is similar to an Einstein black hole potential in four dimensions with 
a cosmological constant, apart from the fact that the Newtonian potential 
is now 1/y/r rather than 1/r! Furthermore, if we were to compare it with 
the six dimensional Einstein black hole (1/r^), we see that gravity in the 
Born-Infeld case is far weaker as one approaches the singularity. 

Keeping this comparison in mind, the flat k = vacua of this theory, 
exist for E = —1 and are non-compact. In order for the region r > rh to 
be spacelike, we need a > and thus a^ < (similar to a planar AdS black 
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Figure 6: The horizon positions r_ and r+ as a function of the black hole 
mass /x for the Born-Infeld case E = +1 and for a < 0, where the vacua are 
dSii. For ^ = ^AT we obtain the degenerate horizon case. The outer horizon 
r_|_ is also the horizon that we obtain for the case E = +1 and for q > 0, 
where the vacua are AdS^^. 

hole). The horizon position is simply r^ = (4a/i)^'^. We will see more of 
the flat vacua in the last subsection. 

For the case E = —1, the potential has only one root. In order for the 
region r > r/^ to be spacelike, we need to have a > 0. Then both cases 
K, = ±1 are possible. For the AdS^ (k = — 1) vacua, the horizon distance is 
bounded from below r^ > 2a, while for the dS^ {k = +1) vacua, the horizon 
is not bounded. The plot of the horizon positions as a function of /x is given 
in Fig. 5. The situation is similar to an AdS Schwarzschild black hole with 
planar or spherical horizon. 

In the E = +1, the potential has either none or two roots. In order for 
the potential to acquire two roots the black hole mass should be in the range 
fi > fiN = 4(2/5)^'^|ap^. Then we can distinguish two possible cases, one 
with a > and n = —1, and a second with a < and k = +1. For the 
AdS4 {k = —1) vacua, the spacelike region is for r > r^ and thus the space 
is non-compact. The horizon is bounded as 2a/5 < r^ < 2a. For the dS^ 
(k = +1) vacua, the spacelike region is for r^ < r < r^^- and thus the space 
is compact. The outer horizon is bounded as 2a/5 < r^ < 2a and the inner 
as < r^ < 2a/5. When lu, = n^ the two roots become degenerate and 
VJ^_^ = 0. This solution will be discussed in the following section. The plot 
of the horizon positions as a function of // is given in Fig. 6. 
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4.2.3 The Born-Infeld Nariai vacuum 

The degenerate case of the Born-Infeld vacuum for k = +1, E = +1 and 
a < deserves special consideration. This happens when fi -^ ^^ and so 
r_ ^ r-i- ^ Y^2|a|/5. The internal space then appears to collapse, however 
one should look as for the Nariai metric in a different coordinate system, as 
we saw in Sec. 4. 1.3. Repeating the same procedure as before and expressing 
the parameters /i and |a| as functions of r_|- and ^ 



r[ 



3 



e(i-e')' , , rii-^/^ 



^ 2 (l-el/2)(l-^5/2) ' l«l 2 1-^1/2' ^^^^ 

we find that the Nariai limit ^ ^ 1 gives a non-singular limit for the modified 
potential / -^ -2{1—B?). Therefore, we obtain a non-vanishing internal space 



ds^ 



2\a\ 



1 ,^ ^os ,„9 dR 



2 



a - R')dQ' + y- — + h^.dx^'dx^ 



(40) 



Furthermore, the deficit angle from (21) can be expressed as a function 
of the modified potential as 

(i^ = ii^ = \\dnnC = \c . (41) 

4.2.4 Flat vacua 

Let us now study the flat (k = 0) vacua that exist in the a^ 7^ case. We 
have found previously a flat vacuum in the special BI case. More generally, 
we know from (15) that flat vacua can exist only for e = —1. The horizon 
position is at r\ = —[ija^ and we can see from the potential asymptotics 
that the region r > r^ will be spacelike in two cases: (i) a < and a^ < 
and (ii) a > and — l/(4a) < a^ < 0. Thus, for both cases we should have 
^ > 0. In case (i) there is a branch cut singularity r^ given in (12) hidden 
always beyond the horizon. In case (ii) there is not such r^. 

In all these cases, the flat vacua have non-compact internal spaces and 
are warped. The latter property is in contrast to the flat vacuum that exists 
when a} = 0. Thus, we see that the introduction of a bulk cosmological 
constant can be compensated by the black hole mass ^i and warping of the 
four-dimensional space, giving again vacua of zero four-dimensional curva- 
ture. As noted before, to compactify such flat vacua, a gauge field flux has 
to be added along the lines of Appendix A. 
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5 Self-properties of the solutions 

Let us now discuss the physical consequences of the above solutions. In par- 
ticular, we wish to see whether we can obtain codiniension-2 braneworlds 
exhibiting self- accelerating or self-tuning behaviour. For the former, we seek 
(iS'4 models where the acceleration is mainly of geometrical origin rather 
than due to the vacuum energy or tension T of the branes. For the latter, 
we seek braneworld solutions (de Sitter or flat) where the effective Hubble 
parameter of the brane is not related to the brane vacuum energy T and the 
tension can be arbitrary large even for small observable Hubble curvature. 
Here it is important to emphasize, that the Gauss-Bonnet term plays the 
role of an induced gravity term on the brane [23] according to (22). There- 
fore even if the internal space is not compact we expect that gravity will 
be four dimensional up to some distance scale. In fact, we would expect, 
and this turns out to be true, that self-acceleration will hold for warped 
but uncompactified codimension-2 setups as for the case of codimension-1. 
Secondly, taking into account that the effective curvature of the brane is 
given by 

i.e. as the sum of a geometrical term and a brane tension term, the self- 
accelerating and the self-tuning cases are, in some respect, in contrast to 
each other. 

Let us first study self- acceleration. The key relation is (42). For this 
case, one should firstly have a non-zero positive geometrical contribution 
to the curvature, i.e. a < coming from the Gauss-Bonnet term in the 
action (1). Secondly, this contribution should be dominant in comparison 
to the brane tension contribution. Note that in order to be in accord with 
phenomenology from supernovae data [1] and since Hq ~ lO^^^ey, the 
Gauss-Bonnet coupling should be roughly of the order a ~ 10^'^^Mp^ . In 
this case, the Gauss-Bonnet term in the action is dominant in comparison 
to the Einstein-Hilbert term. Furthermore, for this limit, even the compact 
models will have enormous volumes of cosmological size, and therefore will 
be comparable to certain aspects of infinite volume ones (see Appendix C 
for the calculation of the volumes of the various models). 

It is useful to define a self-acceleration index as the fraction of the geo- 
metrical acceleration to the curvature 

Vm)^^ (43) 

H^ 2\a\ ■ ^ ^ 
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When < s < 1 we have T > assuming G4 > 0. On the other hand, if 
s > 1, we have T < with the hniit \T\ -^ 3/(167rG4|Q|) giving s -^ 00. 
We can argue that self-acceleration exists whenever 1/2 ^ s < 3/2, with the 
strict self-accelerating limit being s ^ 1 for T = 0. 

The index (43) depends on the specific bulk solutions that we have. 
Going through the dS^ vacua that we discussed in the previous section, we 
should restrict ourselves to the ones with a < 0. Then the self-acceleration 
index for the various cases is listed as following 

1. For a^ = 0, e = +1, a < 0, we have two horizons and the index for 
the two branes varies as < s_ < 1/6 and 1/6 < s+ < 1/2. Thus, for 
this case (Fig.l) we have no strict self-accelerating limit. The Nariai 
limit case gives s = 1/6. 

2. For a^ = 0, e = —1, a < 0, we have one horizon and the index is s > 1 
(Fig. 2). Here we have a strict self-accelerating limit when ^ ^ fig, 
however, at /U = /Ug the model is singular. 

3. For a^ = 0, e = +1, ^ = 0, a < 0, we have the compact Z2-symmetric 
model with s = 1/2. 

4. For a^ 7^ 0, e = ±1, ^ = 0, a < 0, we have the compact Z2-symmetric 
model (Fig. 4). The index for e = — liss>l and for e = +1 varies as 
< s < 1. The strict self-accelerating limit happens for the BI case 
and can be approached by either branch e = 1 and e = — 1. One can 
argue that the BI limit is natural in the sense that the theory becomes 
more symmetric at that point. 

5. For the BI case with E = +1, ^u 7^ 0, q < 0, we have two horizons 
(Fig. 6) and the index for the two branes varies as < s_ < 1/5 and 
1/5 < s+ < 1. The strict self-accelerating limit happens for the ^ = 
case as mentioned above. 

From the above we see that the two physically interesting cases where 
the late acceleration of our Universe can be of purely geometrical origin, 
are cases 2 and 4. In case 2 the strict s = 1 limit can only be reached 
asymptotically in order not to hit the singularity at fi = fig- 

Actually, the first two cases present some similarities and one important 
difference with the codimension-1 DGP [8], which is worth noting. In fact, 
we see that case 2 corresponds to the self-accelerating branch of DGP [14] 
and case 1 to that of the normal branch with ^5*4 branes [11]. A first sim- 
ilarity is the fact that we obtain branching, which now is in-between the 
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Gauss-Bonnet branch and the Einstein branch. Secondly, for case 1, as for 
the normal branch with dS'4 branes, we expect a normalisable spin-2 gravi- 
ton since the relevant volume element is finite. It is an intriguing difference, 
however, that the codimension-1 unstable or self-accelerating branch corre- 
sponds to the Einstein branch here and not the Gauss-Bonnet one. 

In case 4, we note that the limit s = 1 is regular and we see that around 
a^ = l/(2|a|), self-acceleration is possible for both branches. Here, the 
similarity with codimension-1 resides that in that the case of T = tension 
we have enhanced symmetry [10,11], as one also expects for the Born-Infeld 
case in codimension-2. 

Let us now discuss the cases which have to do with self-tuning. The 
self-tuning idea aims to allow for vacua where their effective cosmologi- 
cal constant is independent from (or at least not strongly dependent on) 
the vacuum energy of the brane, without fine-tuning. In the toy-model 
we are presenting here, vacuum energy is represented by the brane tension 
T± whereas the effective cosmological constant is represented by kH±. For 
self-tuning to operate, there should be enough integration constants allowed 
to be adjusted when we vary one brane tension, with the crucial demand 
of keeping the curvature of that brane constant. In the solutions that we 
have discussed, the free parameters are the angular coordinate range c and 
the black hole mass /U. On the other hand, the action parameters are the 
Gauss-Bonnet coupling a, the bulk cosmological constant a^ and the brane 
tensions T±. If we vary for instance T+ we do not wish a, a^, T„, H^ to 
change, but only possibly c and /i. We will be obviously interested in the 
(iS'4 and the flat vacua to test if such a self-tuning can work. 

For the (iS'4 vacua, the Hubble parameter on the brane depends on the 
black hole mass n and the action parameters a and a^. Since we require 
that the action parameters are not tuned, if we fix the curvature H^ this is 
equivalent to fixing the black hole mass fi. In the cases where we have more 
than one brane with different Hubble parameters H±, both of the latter are 
given functions of q, a^ and jj,. Thus, fixing /j, from the curvature of the one 
brane H^ , fixes also the curvature of the second brane H^ . Then from the 
relation 

,.„_,„(_L,,H|)^l=^r^ ,44, 

we see that, if we change T~^, we can keep the curvature H^ constant by 
changing c. But since H- is fixed, we have to change also T~. This results 
to an interbrane fine-tuning. Hence, the only way to obtain self-tuning is 
when the second brane is absent, or in Z2-symmetric models where only one 
of the above relations survive. 
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These selftuning cases, however, are not aU theoreticaUy satisfying if 
we look at the orders of magnitude of the various dimensionful quantities in 
(42). Although we have no idea of the brane tension value today, it would be 
logical that, if no accurate cancellations happen during phase transitions in 
the history of the Universe, its natural order of magnitude is T ~ Mp^. On 
the other hand we know that the the present value of the Hubble constant 
is approximately Hq ~ 10~^^ Mpi. In order that (42) can then be satisfied, 
today's curvature should be negligible in comparison with |a|~^. In other 
words, we should find models which allow for |a|iJ^ <C 1 and approximately 

5^r.J-. (45) 

3 2a 

The above relation should not be viewed as a fine tuning between the action 
parameters a and T, since G4 scales as ^ (23) and can vary in time when 
the vacuum energy changes, in accordance with the remarks of [35]. The 
constraints on the time variation of Newton's constant come rather late 
in the history of the Universe, certainly after the QCD phase transition. 
Therefore, the relation (45) is providing the angle deficit/excess parameter 
P, given a Gauss-Bonnet coupling a, which needs not to be unnatural as 
in the self-accelerating case. Furthermore, the relation (45) is exact for the 
case of flat «; = vacua. 

From the models that we have discussed, let us see the cases that give 
rise to such theoretically viable self-tuning 

• For a^ = 0, e = — 1, and for both a > and a < 0, we have the 
non-compact vacua with |a;|-?/^ <^ 1 for ^|a|~^'^ ^ 1 (see Fig. 2). 

• For a^ / 0, e = — 1, ;U = 0, a < 0, we have Z2-symmetric vacua with 
lalH"^ < 1 for |a|a2 < 1 (see Fig.4). 

• For a^ / 0, e = —1, ^ = 0, a > 0, we have non-compact vacua with 
aH"^ < 1 for aa'^ < 1 (see Fig.4). 

• For the BI case with E = —1, a > 0, we have non-compact with 
aH'^ < 1 for ^a"3/2 ^ i (ggg Fig.5). 

• All the flat non-compact vacua, for both a > and a < are self- 
tuning in the exact sense. 

Here we should note once more that the above (exact or approximate) 
self-tuning solutions do not constitute a resolution, but rather a potential 
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amelioration, of the cosmological constant problem. We have found vacua 
where the brane curvature is independent of the brane vacuum energy and 
furthermore the these two energy scales can be well separated. However, 
the important question, for which the self-tuning mechanism cannot give a 
simple answer, is if the dynamical variation of the vacuum energy results 
in remaining in these vacua. In some sense, the term self-tuning itself is 
deceiving since the existence of these vacua does not guarantee that there 
are some attractor dynamics which tune the system towards them. Never- 
theless, the very existence of them is important, since at least at the level of 
non-dynamical solutions, dissociating the brane curvature from its vacuum 
energy is a rather unique situation. Moreover, we saw that most of these 
vacua are non-compact, but, as we will explain in the last section, gravity is 
quasi-localised and thus effectively four dimensional up to some cross-over 
scale. 

It is evident that the vacua with the above self-tuning properties are 
different from the ones with the self-accelerating properties corresponding 
in particular to a totally different bare parameter a. This is true as long 
as we want self-acceleration at a very low energy scale. If, for example, 
we want to explain inflation theory geometrically, as in [36], then the self- 
accelerating (inflationary) vacua can be deformed to self-tuning ones by 
letting the integration constant ^ run. A nice toy model scenario for this 
is pictured in Fig. 2 if we assume that fi is an increasing function of proper 
cosmological time. One starts at the big-bang singularity which corresponds 
to a 6 dimensional bulk naked singularity at fi = fig. At early time we have 
H^ lis and the inflationary expansion is, Hf^r ~ 1/(2|q|). At late time the 
late expansion is related to the large mass of the soliton whereas the vacuum 
energy is of the order of l/(G'4|a|) which is much larger than the current 
Universe curvature H^Mp^. 

Before we close off this section, it is worth pointing out that no self- 
tuning cases are possible for the Gauss-Bonnet branch. Furthermore, the 
only finite volume self-tuning case that we encountered is possible for a de 
Sitter braneworld. 

6 The Gauss-Bonnet instantons 

In our previous analysis, we have found several brane world vacua of finite 
volume. These solutions are of special interest because they can give rise to 
gravitational instantons. We saw how they all have the feature that the four 
dimensional space is ^5*4. Wick rotating to a Euclidean metric h^^ as in (5) 
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with K = +1 and keeping the Euchdean internal space, all the coordinates 
are spacelike and the solutions describe a Riemannian manifold with conical 
singularities. Particularly important are the instantons which are regular, 
in other worlds those having no conical singularities (see also [37,38]). We 
can divide these vacua into two classes according to their topology. 

The first class of instantons are the ones for vanishing black hole mass 
/i = 0, where there is no singularity at r = and the space is Z2-symmetric. 
These instantons are given by 



ds'={l- ^- de' + -, —^ + r'h^^dx^'dx'' , (46) 



with r^ = — 2a/(l + e\/l + Aaa'^) and are topologically S^. In order for 
these instantons to be regular, one should have f3 = 1, which from (21) gives 
c = rh for the periodicity of the ^-coordinate. These instantons have an 
Einstein theory limit only in the case e = — 1 and 1 + 4aa^ 7^ 0. 

The second class of instantons are the ones of the Nariai limits of vacua 
with non-vanishing black hole mass found in Sec. 4. 1.3 and Sec. 4. 2. 3. These 
instantons are given by 



ds = ri-i^{a, a ) 



dR^ 

C(l - i?2) "^ '"^-^^ 



C(l - R^)de' + — -J- + h^udx^'dx' 



(47) 



where (^ is a numerical factor and r^ the degenerate horizon. These instan- 
tons, in contrast with the ones noted before, are topologically S*^ x S'^. In 
order for these to be regular, one should have (3 = 1. The deficit angle for 
the instantons that we have found is given by /3 = C,C, thus, regularity of 
the instantons imposes C = 1/C,. Let us also note here that these S*^ x S'^ 
instantons for the case a? 7^ 0, exist also for the case other that the BI 
limit that we have discussed in Sec. 4. 2. 3. This is why in (47) we allowed for 
a^ dependence of the degenerate horizon. None of these instantons has an 
Einstein theory limit. 

The above instantons have the physical interpretation of describing prob- 
abilities of nucleation processes [39, 40] between two distinct gravitational 
backgrounds. In particular, the probability of nucleation of a pair of Nariai 
black holes from the pure dS^ backgrounds of /i = 0, is given by 

r = r/exp[2A5] , A5 ^ 5^6 - 832^3^ , (48) 

where 77 is the one loop contribution from the quantum quadratic fields and 
556 and Sg2 y. g4 are the values of the action for the two instantons that we 
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presented earlier. It is straightforward to compute the above probabihty of 
nucleation for the two simple cases of a^ = and the BI case a} = l/(4|a|). 
In both cases one obtains 

A5 oc M^a^ > . (49) 

Since these action differences are positive definite, we deduce that the prob- 
ability of both processes is unsuppressed unlike the case of [39]. Therefore, 
the solitonic vacua (^u 7^ 0) are apparently stable with respect to pure de 
Sitter vacua in Lovelock theory. This result certainly deserves further study. 



7 Discussion and Conclusions 

In this paper we studied in some detail codimension-2 braneworlds in the 
context of Lovelock gravity. Using a modified version of Birkhoff 's static- 
ity theorem [24,25], we found all six dimensional solutions describing a de 
Sitter, flat or anti de Sitter braneworld of codimension 2. First thing we 
can observe is that, unlike the codimension- 1 case where the Gauss-Bonnet 
invariant plays a somewhat secondary role, here, in the case of codimension- 
2 braneworlds, it can give rise to self-acceleration and certain self-tuning 
properties which are not present in the Einstein theory. This is largely due 
to the fact that codimension-2 junction conditions induce on the brane an 
Einstein-Hilbert term [22]. It is important to note that this does not mean 
that we have ordinary four-dimensional gravity, rather, as noted in [23], we 
are in a quite similar situation as for the codimension-1 DGP model [8]. By 
this, we mean that from some ultraviolet scale up to some infrared scale 
we expect gravity to "look" four-dimensional. A proof of this statement 
of course requires the full spectrum of gravitational perturbations for the 
solutions that we have found. 

The important relation we came across is (25) which relates the vacuum 
energy or tension T of the brane with the Gauss-Bonnet coupling a and the 
effective cosmological expansion on the brane K^^ which itself is related to 
the characteristics of the bulk solution: its mass, bulk cosmological constant 
(and charge) in particular. The interesting feature we found here is that the 
Gauss-Bonnet coupling can give de Sitter branes without vacuum energy 
on the brane, purely geometrically. In this sense such de Sitter solutions 
are self-accelerating. In order to explain the small cosmological constant we 
then have to fine-tune as usual a ~ H'^ . 

The second important point is that Gauss-Bonnet coupling a and the 
topological parameter /3, which is otherwise unconstrained, dictate the cross- 
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over scale between the four-dimensional Planck scale G4 and the six-dimensional 
one Gq (23). Here, unlike the five dimensional DGP model, the cross over 
scale is not necessarily tied up to the self-acceleration scale. Indeed, the 
more (3 is close to 1 (23), the more we can dissociate these scales. In a nut- 
shell, all depends on the gravity phenomenology beyond the cross-over scale 
and is dictated by the full graviton propagator on the brane. As we empha- 
sized, the appearance of the induced Einstein tensor in the codimension-2 
junction conditions is not a proof of an effective four-dimensional gravity or 
not. It is a definite sign, however, that a relevant scale will appear in the 
boundary conditions for the gravitational perturbations. Although we have 
not undertook the precise perturbation calculation here, we will comment 
on some of its characteristics later on. 

The most conservative approach is, as in DGP [14], to introduce a hier- 
archy of scales between Gq and G4 so that the cross-over scale is 

r'. = ^- (50) 

(-T4 

Assuming that re ~ H^ , this gives a very low higher dimensional Planck 
scale Me = Gq^^"^ ~ W'^^Mpi which in turn dictates that 1 - /3 ~ 0(1). 
But, as mentioned above, one can have re <C H^ by having 1 — /3 ~ 0. For 
the self-tuning vacua, on the other hand, the natural value of the Gauss- 
Bonnet coupling is a ~ ^pi which for the same cross-over scale (50) and 
Tc ~ Hq gives a huge angle excess 1 — /3 ~ lO^^'^. Obviously, in this 
case a reasonable hierarchy between Hq and the cross-over scale r^, cannot 
reduce significantly the latter huge excess angle. We emphasize, however, 
that the above orders of magnitude, should be viewed with caution since the 
definition of the cross-over scale (50) should be done in a true cosmological 
setting with complete knowledge of the modified FRW equations. 

Another important point concerns the fact that the self-accelerating 
branch of DGP seems to be embedded in the usual Einstein branch of Love- 
lock theory and not in the Gauss-Bonnet branch. Indeed, given the Gauss- 
Bonnet term in the action and no bare cosmological constant (a^ = 0) we 
find two distinct solutions for e = ±1 (7): in the Gauss-Bonnet branch, 
e = +1, we can have a finite volume soliton solution since then the Gauss- 
Bonnet coupling plays the role of a cosmological constant (Fig.l) (situation 
akin to a Schwarzschild de Sitter black hole [31]). Then, we have a finite 
volume element and we expect a normalisable zero mode graviton. We get 
no self-acceleration. For the Einstein branch however, e = —1, we have a 
single brane, infinite volume element and self acceleration with as small ten- 
sion as we want (Fig. 2). The bulk solution corresponds to a Gauss-Bonnet 
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corrected Schwarzschild black hole. This difference may be interesting in re- 
spect to linear (in)stability and the presence of ghosts of such backgrounds, 
not only in the scalar, but also in the spin-2 sector. 

Furthermore, in our solutions we see at least two different length scales 
emerging where we can probe differing four-dimensional gravity. That of 
the volume element, which if finite, would mean that beyond that scale we 
expect ordinary four-dimensional gravity. Secondly, that of the cross-over 
scale where up to that distance gravity seems four-dimensional, as in the 
case of the codimension-1 DGP model. Again we emphasize that the cross 
over scale and self acceleration scale are here supplemented by an extra 
topological parameter /? which in a sense tells us how far we are from a 
Kaluza-Klein setup for the Killing direction dg. Indeed, we see that for 
/? ^ our codimension-2 space is very much like a lightning rod [41], where 
we seemingly have a cascade from a six-dimensional to a five-dimensional 
and then a four-dimensional setup (reminiscent of the recent proposal [42] ) . 

What is missing from our analysis is the linear perturbation of these 
solutions which will tell us of the stability and the precise gravitational 
spectrum. This is not a trivial task, for our metric is not conformally flat. 
Therefore, apart from the usual complications of black hole perturbations, 
one has to add the fact that Lovelock perturbations are going to be genuinely 
different from the usual ones for Einstein theory. The reason is simple: 
the background Weyl tensor appears in the Lovelock field equations and 
thus, extra tensor pieces are bound to appear for tensorial perturbations. 
Recent work in fact from string theory [43], in relation to the AdS /CFT 
correspondence, indicates that there is yet another scale appearing in the 
bulk perturbations of (6). This extra scale comes from the extra tensorial 
pieces "to be added" to the usual perturbation operator, that dictate that 
gravity waves in a planar (k = 0) black hole background (6) can evolve in 
a differing background geometry than that of (6). This means in particular 
that the light-cone of the wave fronts can break causality. We think that 
this is an intriguing and extremely important issue, which we hope will be 
undertaken soon. 

In addition, the cosmology of codimension-2 branes is very poorly un- 
derstood, especially from the point of view of self-tuning and the issue of 
vacuum selection. It has been known that in Einstein gravity the intro- 
duction of matter other than tension on codimension-2 branes introduces 
singularities'^ worse than conical [45] (see, however, the work of [46] for in- 



^This is a well known fact in four dimensional gravity, called the Israel paradox con- 
cerning self-gravitating cosmic string metrics (see for example [44]). 
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tersecting brane cosmology in six dimensions). Although, the appearance 
of singularities is natural in defects of codimension higher than one [47,48], 
one has in practice to regularise the brane [49, 50] in order to consider some 
cosmological fluid on them [51,52]. Although one expects that brane singu- 
larities persist in the Lovelock theory, the fact that there exists an induced 
Einstein equation on the brane, allows for the cosmology to be studied with- 
out the need of explicit regularisation. We plan to address the question of 
cosmology of the present Lovelock models in a different publication [53] . 
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Appendix A: A more general action 

In this Appendix we will consider a generalization of the model that we 
discussed in the main text to a D-dimensional bulk spacetime, where also a 
gauge field Fmn is coupled to gravity. The action of the system reads 



D, 



S= d^'x V=^ 



{R + aCcB) - A 



IQttGd 



\jd^x,/^gF\ (A.l) 



where Gd the D-dimensional Newton's constant. 

It is straightforward to write down the Einstein equations of motion for 
the above action. They read 

Gmn — O.HMN = SttGdTmn , (A-2) 

with the energy momentum tensor 

Tmn = —J^guN + FmkFn ~ J^ 9MN ■ (A-3) 

Furthermore, the gauge field equation reads 

dMiV^F^'') = 0. (A.4) 
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The static spherically symmetric solutions of the above equations of mo- 
tion with {D — 2)-dimensional space of maximal symmetry are locally iso- 
metric to 

flv 
ds^ = -V{r)dt^ + —— + r^K^dx^'dx" . (A.5) 

V[r) 

The explicit solution of the equations of motion then gives 

q 

47rr 



^ = TZ^dt A dr , (A.6) 



and for the potential, 
Vir) = .+ ^ 



l + eWl + 4a(a^ '^ ^' 



(A.7) 



(A.8) 



^D-l ^2(D-2) 

and the parameters appearing in the action have been rescaled to, 
a = {D - 3){D - 4)d , WttGdA = {D - 1){D - 2)a^ (for dSo) 

The integration constants are, 

p2 _ Gp q^ _ IGttGdM 

^ 27r(Z)-2)(Z)-3)a ' ^ (D - 2)S« ' ^ ' 

where q is the charge, M is the AD or ADM mass of the solution and S^ is 
the area of the unit {D — 2) maximally symmetric subspace. Finally, e = ±1 
which gives rise to 2 branches of solutions. The convention of the /i sign is 
chosen so that the gravitational mass is always /U > 0. As one can easily 
check by expanding the square root for large distances, the sign flip in front 
of fi is necessary to match the Schwarzschild solution behaviour for positive 
AD mass. On the other hand, the charge term of the above potential is the 
opposite of a charged black hole for e = — 1 . 



Appendix B: The general solution of x^ + Ax + M = 

To determine the positions of the horizons for the a^ = case, we need 
the general solutions of the third order algebraic equation 

x^ + Ax + M = 0, (B.l) 

The solutions of the above equation depend on the parameters A, M and in 
particular on the discriminant of the system 
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Figure 7: The positive real roots of the equation x^ + Ax + M = 0. The 
lower curve gives the root for the ^ > case and the upper one the root(s) 
for the ^ < case. In the interval < M < Mc, there are two roots r_ < r+ 
which become degenerate for M = M^.. 



If Z) > there is only one real root. If Z? < there are three real roots, two 
of which become degenerate where the inequality saturates. A critical value 
of M that will be useful in the subsequent study is 

2 



Mr 



3^/3 



A\ 



(B.3) 



Since the solutions of the above equation will correspond to horizons, we 
also need to look which of these real roots are in addition positive. All 
permissible cases are listed as following 



If ^ > and M < 0, it is D > and we have the root 



rh 



" + ./nl 


1/3 


A 


2 J 




3 



M r- 



"1/3 



(B.4) 



If ^ < and M < —Mc, it is -D > 0, but we have one positive root 
which is given by the same formula as above (B.4). [It has different 
value, however, since A < 0.] 



• li A < and < M < Air, it is D < and we have two roots 



1^+ 
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COS I - COS 
o 



+V3 



sm 





Ml 


\ 
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(B.5) 



(B.6) 
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• If ^ < and —Mc < M < 0, it is I? < but we only have one positive 
root, given by (B.5). 

The results of the above analysis are summarized in Fig. 7. 



Appendix C: Volume calculation 

In this Appendix, we will calculate the volume of the general brane 
world models that we considered in the main text. To do this, we will make 
the assumption that the zero mode graviton wavefunction follows the warp 
factor. In other words, we will assume that the zero mode is emanating from 
the metric 

V{r) 
with the four dimensional metric gfiJ (x) = h^u + H^^, where H^i, is a 



ds^ = V{r)de^ + -^ + r'^gl^J{x)dx^'dx'' , (C.l) 



perturbation around the background hf^^. To find the volume of the space 
we should substitute the metric (C.l) in the action (1) and integrate out the 
internal space. We will define the volume (Vol) to be given by the coefficient 
of the four-dimensional curvature term in the action as^ 



s- (™> 



levrGf 



f d'^x./^^R^ + ... , (C.2) 



where the dots represent higher order in curvature terms. It is straightfor- 
ward to compute that the various quantities of the action (1) for the ansatz 
(C.l) give 



\/^ = r y/^g2 , (C.3) 

i?6 = ^^4 + ... , (C.4) 

^GB = ^(-] Ra + ... (C.5) 



Then the volume is given by the following integral 

'V 



(Vol) = 27rc / dr r^ 



1-F4a 

r 



(C.6) 



"Notice that according to the ansatz (C.l), it is [r] = {d\ = M^^ and thus [x^^] = M'^ . 
For this reason, we obtain an unusual mass dimension for the volume [(Vol)] — M~ . 
With a suitable four-dimensional coordinate rescaling, e.g. x'^ = r^x^, the volume will 
be redefined as (Vol)jv — (Vol)/r^ and will acquire the correct M^^ mass dimension. 
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If 1 + Aaa^ / 0, the above integral gives 



(Vol) 



2ttc 
~3~ 



o 2a / 5/x F 

3 V VI + 4Qa2 r2 



'•+ 



(C.7) 



where F is the hypergeometric function 



iFi 



2 1 7 



Aaefi 



5' 2' 5' (l + 4Qa2)r5 
For the special BI case, 1 + Aaa^ = 0, the integral gives instead 



(Vol) 



27rc 
~3~ 



13 



-r + AaKr 



r+ 



(C.8) 



(C.9) 



It is straightforward to see that the above formulas give finite (and posi- 
tive) volume for the two Nariai cases. To obtain the volume in this case one 
should move to the {R, 0) coordinates and take the limit ^ ^ 1. 
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